Non-LTE Effects of Na I in the Atmosphere of HD209458b 
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ABSTRACT 

The recent announcement that sodium absorption has been observed in the 
atmosphere of HD209458b, the only EGP observed to transit its parent star, 
is the first direct detection of an EGP atmosphere. We explore the possibility 
that neutral sodium is not in local thermodynamic equilibrium (LTE) in the 
outer atmosphere of irradiated EGPs and that the sodium concentration may be 
underestimated by models that make the LTE assumption. Our results indicate 
that it may not be necessary to invoke excessive photoionization, low metallicity, 
or even high altitude clouds to explain the observations. 

Subject headings: planetary systems — radiative transfer 



1. Introduction 

HD209458b is the first extra-solar giant planet (EGP) observed to transit its parent 
star (Charbonneau et al. 2000; Henry et al. 2000) and, consequently, its mass (0.69 Mjupiter) 
and radius (1.35 Rjupiter) are known to a high degree of accuracy. These results leave little 
doubt that HD209458b is a gas giant, and its special orbital inclination has drawn a great 
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deal of attention from observers and theorists over the past year. Recently, an increase in 
the sodium absorption (relative to the continuum) at 5893A was observed with the Hubble 
Space Telescope (HST) during several transits of HD2094598 (Charbonneau et al. 2001). 
The additional sodium absorption is believed to be due to Na D absorption in the EGP's 
atmosphere as stellar light passes through the planetary limb. This observation marks a 
major turning point in the study of EGPs, for now we have direct evidence of an atmosphere 
around HD209458b and a measurement of one chemical species (Na). According to Char- 
bonneau et al. (2001), the HST observations suggest that either the EGP atmosphere has a 
low concentration of neutral atomic Na (due to photoionization, molecular formation, or an 
overall low metaUicity) or that high altitude clouds exist and reduce the amount of stellar 
flux transmitted though the EGP hmb. 

In this letter, we offer an alternate explanation for the observed Na absorption and 
explore the possibility that the Na D feature is altered by nonlocal thermodynamic equi- 
librium (NLTE) effects brought on by the impinging stellar radiation field and insufficient 
coUisional thermalization. If true, NLTE effects would offer a natural explanation for the 
apparently low sodium absorption observed in the HD209458b without the need for excessive 
ionization, a reduced metallicity, or extremely high altitude clouds. Below we will present 
theoretical predictions for NLTE Na D doublet line profiles for the transmitted spectrum of 
an irradiated EGP atmosphere and provide stringent hmits on the non-LTE effects. 



2. Model Construction 

The construction of the model atmospheres presented below follows the procedure out- 
lined in Barman et al. (2001) (here after BHA). In order to make a direct comparison with 
HD209458b, we have adopted the best-fit values for the radius, mass, and orbital separa- 
tion of HD209458b and its parent star pubhshed by Brown et al. (2001) and Mazeh et al. 
(2000) (see table 1). We have also assumed that the total flux received by the planet has 
been uniformly distributed over the planet's dayside. Therefore, the incident flux has been 
weighted by 0.5, whereas in the BHA models, no redistribution was assumed. The intrinsic 
effective temperature of the model (i.e. the effective temperature of the planet in the absence 
of irradiation) is 500K. However, this intrinsic temperature has little effect on the structure 
of the outer atmosphere which is completely determined by the incident flux. With this 
redistribution, our models represent an upper limit to the equilibrium effective temperature 
which is about 1800K. We further assume that the planet has a solar composition with an 
opacity setup identical to the "AMES-cond" models of BHA and AUard et al. (2001). In this 
situation, dust grains form in the atmosphere at locations determined by the chemical equi- 
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librium equations but their opacity contribution is ignored, mimicking a complete removal 
of the grains by efficient gravitational settling. Therefore, the models in this study represent 
cloud free atmospheres. 

We model the flux transmitted through the limb of the planet's atmosphere, by solving 
the spherically symmetric radiative transfer equation (SSRTE) suitably adjusted to account 
for the incident radiation. The incident flux is taken from a separate calculation which 
reproduces the observed spectrum of HD209458. The SSRTE has a significant advantage 
over the plane-parallel solution because a geometry more appropriate (and more accurate) 
for the upper atmosphere is already incorporated. When solving the SSRTE, the atmosphere 
is modeled as a discrete number of concentric shells surrounding the interior (or "core"). The 
solution along characteristic rays that pass through the outer most shells is all that is needed 
to calculate the transmitted intensities and automatically accounts for the curvature of the 
limb. The total transmitted flux is obtained by integrating over the planet limb: 



where r is the perpendicular distance from the planet center to a tangential characteristic 
ray. -Rmax is the planet radius at Tstd = and Tgtd is the optical depth at 1.2fim. Rmin is 
chosen such that for r > Rmin and A = 5880 A, /(r) > 0.1/*(r) where /* is the incident 
stellar intensity. This lower limit ensures that the transmitted intensities are well sampled 
and light from the core region is excluded. Note that both i?max and -Rmin are determined 
by the simulation and depend on the resulting structure of the model atmosphere. The only 
prescribed radius is Rp — r(Tstd = 1), which is equal to the planet radius given in table 1. 
With these definitions, the thickness of the hmb {H — Rmax — Rmin) is roughly 0.07i?p. 



When a gas is assumed to be in LTE, the level populations for each species depend en- 
tirely upon the gas temperature and electron pressure and are given by the Saha-Boltzmann 
distribution. In general, the LTE assumption is a matter of computational convenience and 
is not expected to be valid in most cases, especially in the optically thin regions of an at- 
mosphere. However, it is often assumed that LTE is achieved in very late type stars (even 
Brown dwarfs) despite the fact that this assumption has not been thoroughly tested. De- 
partures from LTE have previously been investigated for Ti I (Hauschildt et al. 1997) and 
CO (Schweitzer et al. 2000) for cool M dwarfs and for CH4 in the Jovian planets (Appleby 
1990). These earher works found that, for these particular species, NLTE effects were small. 




(1) 



2.1. 



Na in non-LTE 
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Due to the close proximity of HD209458b to its parent star, the planet's atmosphere 
is subjected to intense stellar radiation. This inherently nonlocal source of radiation dra- 
matically alters the conditions of the outer atmosphere compared to an isolated EGP. In 
addition, the cool atmospheres of EGPs are dominated by strong opacity sources like H2O 
and CH4 and, therefore, have intensities that differ greatly from those of a blackbody. When 
combined with the relatively low pressures of the outer atmosphere, these conditions very 
likely will lead to departures from LTE which may play an important role in determining 
the atmospheric structure and resulting spectrum. If LTE is to be achieved in EGP atmo- 
spheres, then the coUisional rates must be large enough to compensate for the deviations of 
the radiative rates from their LTE values. Unfortunately, a hmitation of all NLTE models 
is a lack of well determined collisional cross-sections for interactions with important species 
like H2, He, and H. However, it is unlikely, given the low thermal velocities, that these parti- 
cles could restore LTE completely. In fact, the effects of collisions with hydrogen on the Na 
D profile in M dwarfs with chromospheres arc fairly small and electron collisions dominate 
despite the fact that Nn/Ne ~ 10^ (Andrctta ct al. 1997). 

The Na I model atom used in this work includes 53 levels, and 142 primary transitions 
(all bound-bound transitions with \og{gf) > —3.0) were included in the solution of the sta- 
tistical equilibrium equations. For details of the method used to solve the rate equations, see 
Hauschildt & Baron (1999). As in Hauschildt et al. (1997) for Ti I, electron-impact bound- 
free collisional rates are approximated by the formula of Drawin (1961) and bound-bound 
collisional rates are based on the semi-empirical formula of Allen (1973) with permitted 
transitions determined from Van Regemorter's (1962) formula. For the ground-state pho- 
toionization cross sections, we have used the Opacity Project data of Bautista et al. (1998). 
We have also assumed complete redistribution. 

We have constructed four model atmospheres and resulting transmitted spectra for 
HD209458b. The first irradiated model (A) only includes collisions with electrons, where 
the source of free electrons is primarily the ionization of potassium. Even in an irradiated 
EGP atmosphere, the number density of electrons is very small {Ne/Ni{2 ~ 10~®) and the 
temperatures are too low for electronic collisions to be important. Model A, therefore, can 
be considered a lower limit for the collisional rates. The second irradiated model (B) has 
the same parameters as model A, except that coUisions with H2 are also included. In order 
to place a secure upper limit on the effects of collisions with H2, we have treated H2 as if it 
had the same rate coefficients as an electron. Assuming identical cross-sections implies that 
the rate coefficients scale with (/i)~2, where /i is the reduced mass of the collision partners. 
Therefore, we are overestimating the H2 collisional rates by more than an order of magnitude. 
The remaining two models (C and D) represent non-irradiated atmospheres with the same 
effective temperature as the irradiated models (~ 1800K) and include the same lower (model 
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C) and upper (model D) limits for the collision rates. Each of these model atmospheres was 
produced from a self-consistent solution of the SSRTE, chemical equilibrium equations, and 
the NLTE rate equations. We did not prescribe the temperature structure or mixing ratios 
for any of the species. 

3. Results 

The temperature-pressure profiles for models A and C are shown in figure 1. Note that, 
the profile for model B is identical to that of model A and model D has the same profile 
as model C. The gas pressures in the limb are very low, ranging from 0.004mb to 0.2mb. 
As is to be expected, the non-irradiated models look nothing like the irradiated models 
which have temperatures nearly lOOOK higher in the outer atmosphere (Barman ct al. 2001; 
Goukenleuque et al. 2000; Seager & Sasselov 1998). There were no significant changes to the 
temperature-pressure profiles in the four NLTE models compared to the LTE structures. 

Departures of the level populations from LTE, for a particular species, are usually 
described by the departure coefficients, bi = n*/ni, where n* is the NLTE population density 
for level i and n-j is the LTE value (Mihalas 1970). In figure 2, we show the 6i's for the 
ground state (3s) and the first excited states (3pi, Spa) of the neutral sodium atom for 
several different physical conditions within an EGP atmosphere. The largest departures 
are seen in model A where the 3s (ground state) and 3p levels are both underpopulated 
by many orders of magnitude, especially in the limb region. The main reason for such 
large departures is that we have only included collisions with electrons and, in such a cool 
atmosphere, the number density of electrons is 8 orders of magnitude below that of 
the dominant species, H2. Consequently, there are essentially no collisions in model A to 
thermalize the level populations in the Na atom thus allowing the radiative rates to dominate 
and drive the system out of LTE. Also, since the upper atmosphere is dominated by a large 
external radiation source (which happens to peak near the Na D doublet, ~ 5000A), the 
mean intensity of the line is much larger than the thermal source function which implies a 
strong decoupling of the radiation field from the local conditions. Furthermore, the ratio 
of the line source function to a blackbody is roughly given by the ratio of the departure 
coefficients of the upper (3p) and lower (3s) levels for the transition (Bruls et al. 1992). 
From figure 2, we see that the line source function is far from a blackbody for the majority 
of the upper atmosphere (since b^p/h^s ^ 1)- 

An obvious question to ask is whether the system will return to LTE if collisions with 
the dominant species (H2) are included. Model B, which has the same parameters and 
temperature profile as model A, includes collisions with both electrons and H2 (but assuming 
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that H2 has the same rate coefficients as an electron). In effect, we have increased the 
importance of electronic collisions by more than 8 orders of magnitude. As is expected, the 
departure coefficients are closer to one, but the level populations in the limb are still far from 
the LTE values. However, since h^p/h^s ~ 1 for most of the atmosphere, NLTE effects will 
only be important in model B for Tgtd < 10~^. The real situation is likely to be somewhere 
between models A and B. 

We also show departure coefficients for non-irradiated models (C and D) with T^q equal 
to the equilibrium effective temperature of the irradiated atmospheres (1800K) for both 
collisional rate limits. Despite the absence of any extrinsic radiation and much lower tem- 
peratures in the upper layers, the non-irradiated atmospheres still have departures from LTE, 
though generally less significant than those found in the irradiated atmospheres. Model C 
(the near coUision free limit) , has departures from LTE similar to those of the coUision dom- 
inated irradiated model (B) for the ground state but with an over populated 3p level. With 
increased collisional rates, Na I has nearly returned to LTE for most of the atmosphere in 
model D with departures still present in the top most layers whereas in the irradiated case 
both levels were greatly underpopulated even in the collision dominated model (B). 

The effects on the Na D fine profiles are quite dramatic for model A (see figure 3). In 
this case, the lack of thcrmalization reduces the line transfer to nearly a pure scattering 
case. As a result, the doublet appears completely in emission. However, in model B, as the 
collisional rates arc increased, the line wings return to their LTE shape while the line cores 
arc reversed. The data analysis of Charbonncau ct al. (2001) does not directly reveal the 
Na line profile produced by the planet atmosphere. Instead, their work shows that a deeper 
transit is observed in the Na band implying additional Na absorption by the planet limb. 
Given the current sensitivity of the observations, a core reversal feature like the one shown in 
our model B would be buried in the noise and manifest itself as simply a smaller equivalent 
width. The Na absorption in model B would result in a transit deeper than in the continuum 
bands but not as deep as implied by the LTE model. The fact that Na D absorption (and 
not emission) has been observed in the transmitted spectrum of HD209458b rules out our 
model A indicating that some thcrmalization does occur in the limb region. However, it is 
unlikely that the collisional rates are as large as those in our model B suggesting that the 
equivalent width of the Na D doublet will be substantially reduced by NLTE effects. 

The reduced equivalent width predicted by our model is not due to photoionization of 
Na. In the majority of the limb, only 3% is ionized and the neutral Na concentration is 
nearly constant with iVNa/A'H2 ~ 10^^'^. Na is only significantly ionized 5%) at the very 
top of the atmosphere where Pgas < lyubar. The shallow ionization depth of Na is due to the 
strong UV opacity provided by metals (e.g. atomic Mg, Al, Ca, Fe, and Ni) which effectively 
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shield Na from the incident ionizing photons. The ionization predicted by our models (even 
3%) is far greater than what is obtained from an LTE calculation but does not significantly 
affect the line profile. However, if the planet's atmosphere is substantially cooler than in 
our model, then additional condensation and settling could further deplete the atmosphere 
of metals and allow greater ionization of Na to occur. We also find that only a very small 
amount of Na is in molecular form and that neutral Na is nearly 3 orders of magnitude more 
abundant than NaCl (the most abundant Na bearing molecule) . Furthermore, condensation 
of Na via NaCl grains is unhkely. The reduced equivalent width and central core emission is 
purely a radiative transfer effect. 

4. Conclusions 

The recent HST observations provided the first direct measurement of the conditions 
inside the atmosphere of HD209458b. However, even under the best circumstances, determin- 
ing the concentration of any species based solely upon one absorption feature is problematic, 
especially if this feature forms in the upper regions of an atmosphere where pressures are low 
and NLTE effects are greatest. Our models clearly show that Na is far from being in LTE 
in the upper atmosphere of HD209458b and the observed Na absorption can be explained 
with a solar metallicity atmosphere which is cloud free or has only very low lying clouds. It 
is likely that other important species (e.g. CO and CH4) are in NLTE, and we plan to test 
the LTE assumption for a wide variety of atomic and molecular species in a future work. 
Hopefully the Na D doublets and other alkali metal lines will be useful diagnostics in the 
study of EGP atmospheres. However, only with detailed NLTE calculations including well 
determined coUisional rates will we have a chance at constraining the physical conditions in 
the atmosphere of HD209458b. 

We would like to thank Charbonneau et al. for making a copy of their work available 
before publication. This research was supported by NASA ATP and LTSA grants to the Uni- 
versity of Georgia and to Witchita State and in part by the Pole Scientifique de Modehsation 
Numerique at ENS- Lyon. 
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Table 1. Model Parameters 





HD209458b 


HD209458 


radius 


1.35 R Jupiter 


1.146 Ro 


mass 


0.63 Mjupiter 


1.05 Mo 


TcfT 


1800 K 


6000 K 


[Fe/H] 


0.00 


0.00 


semi-major axis 


0.045 


AU 
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Fig. 1. — The temperature-pressure profiles for the irradiated (sohd fine) and non-irradiated 
(dashed fine) models. The unshaded area indicates the limb region for the irradiated models. 
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Fig. 2. — The departure coefficients for neutral sodium in tlie atmospliere of HD209458b 
as a function of tlie radial optical depth at 1.2|Lim(rstd) and pressure (top axis). Model 
A demonstrates departures for an irradiated atmosphere including only collisions with free 
electrons. Model B is the same as A, but includes collisions with H2 assuming that H2 has 
the same rate coefficients as an electron. Model C shows the departures in a non-irradiated 
atmosphere (where T^q — T^s) with coUisions treated as in Model A. Model D is the same as 
C but with collisions treated as in Model B. Thick hnes refer to the 3s level (ground state) 
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Fig. 3. — The Na D doublet for Model A (top dotted line), Model B (lowest dotted line), 
and when assuming LTE (solid line). The flux has been normalized to one at 5880A. 



